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Abstract. Young massive stars in the central parsec of our Galaxy are best explained by 
star formation within at least one, and possibly two, massive self-gravitating gaseous discs. 
With help of numerical simulations, we here consider whether the observed population of 
young stars could have originated from a large angle collision of two massive gaseous clouds 
at R ~ 1 pc from Sgr A* . In all the simulations performed, the post-collision gas flow forms an 
inner nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the 
observations. Furthermore, the radial stellar mass distribution is always very steep, E* oc R~ 2 , 
again consistent with the data. The 3D velocity structure of the stellar distribution is however 
sensitive to initial conditions (e.g., the impact parameter of the clouds) and gas cooling details. 
In all the cases the amount of gas accreted by our inner boundary condition is large, enough to 
allow Sgr A*to radiate near its Eddington limit during ~ 10 5 years. This suggests that a refined 
model would have physically larger clouds (or a cloud and a disc such as the CND) colliding at 
a distance of a few pc rather than 1 pc as in our simulations. 



1. Introduction 

The central parsec of our Galaxy is host to young, massive "He-I" stars that can be classified 
into two stellar discs that orbit the central supermassive black hole, Sgr A*, outside the inner 
arcsecond. Most of the stars are located in a well-defined thin stellar disc that rotates clockwise 
as seen on the sky BUG}. The rest of the stars can be classified as a second more diffuse disc that 
is more eccentric and rotating counter clockwise @; S| although the statistical significance of this 
feature is disputed by [4(]. The stellar populations are relatively co-eval, formed approximately 
6 ± 1 million years ago 0]. In terms of formation scenarios, the "in-situ" model is currently 
the one best favoured by the Galactic Centre (GC) community; in this picture, stars form 
inside a self-gravitating gaseous disc [e.g., [H; H; [Bj. Theorists had expected gaseous massive 
discs around supermassive black holes to form stars or planets [e.g., 0; 0] long before the 
properties of the young stars in the GC became known. Recently, g have numerically simulated 
the fragmentation process of a geometrically thin gaseous disc of mass ^ 10 4 Solar mass for 
conditions appropriate for our GC (albeit with a rather simple cooling prescription) and found 
a top-heavy mass function for the stars formed there. [9J extended these numerical studies to 
the fragmentation of eccentric accretion disks. 

However, the in-situ model has not so far addressed in detail the origin of the gaseous discs 
themselves. Similar discs are believed to exist in AGN and quasars, and indeed are invoked 
as a means of feeding a supermassive black hole's immense appetite for gaseous fuel [e.g. [T^. 



However, Sgr A*is currently not a member of the AGN club, as its bolometric luminosity is 
around ~ 10~ 9 of its Eddington limit [e.g.,[ll|], implying a quiescent mode of accretion. 

The viscous timescale of a thin, marginally self-gravitating disc around Sgr A*, given a disc 
mass of ~ 10 4 Mq [reasons for this choice of mass see in, e.g., [l2|] is very long compared with the 
age of the stellar systems (approx. 10 9 yrs compared to 10 6 yrs). In fact, the gaseous discs would 
have to evolve even faster as the two stellar systems are largely co-eval. The requisite gaseous 
discs could therefore not have been assembled by viscous transport of angular momentum. 

There is also no strong evidence for other similar star formation events in the region within 
the last ~ 10 8 years. Therefore, the one-off star formation event appears to be best explained by 
a one-off deposition of gas within the central parsec. There are several ways in which this could 
have happened, e.g. a Giant Molecular Cloud (GMC) with a sub-parsec impact parameter (in 
relation to Sgr A*) could have self-collided and become partially bound to the central parsec 
[e.g., [El]. Alternatively, a GMC could have struck the Circumnuclear Disc (CND) located a few 
pc away from Sgr A*, and then created gas streams that settled into the central parsec. 

In this article we explore such a one-off collision event in a very simple setup. We allow two 
massive, uniform and spherical clouds on significantly different orbits to collide with each other 
one parsec away from Sgr A*, and follow the gas dynamics in some detail. We find that the 
collision forms streams of gas with varying angular momentum, both in magnitude and direction. 
Parts of these streams collide and coalesce to form a disc in the inner region. As the gas cools, 
it becomes self-gravitating, and stars are born. The resulting distribution of stellar orbits is 
compared with the observational data throughout the article, and particularly in ^3 We find 
that our model is reasonably successful in explaining the population of young He-I stars in the 
Galactic Centre. 

2. Numerical methods 

The numerical approach and the code used in this article is the same as in (§[ with only slight 
modifications. Namely, we employ GADGET-2, a smoothed particle hydrodynamic (SPH)/N- 
body code. Gas cools according to du/dt = —u/t coo \(R), where cooling time depends on radius 
as 

^cool(-R) = (3tdyn(R) (1) 

where t dyn = and O = (GM bh /i? 3 ) 1/2 , the Kepi erian angular velocity, where (3 is a 

dimensionless number. Gas moves in the gravitational potential of Sgr A*, modelled as a 
motionless point mass with Mbh = 3.5 x 10 6 Mq at the origin of the coordinate system, and a 
much older relaxed isotropic stellar cusp potential given by the mass profile derived by [2j. 

To model star formation, we introduce sink particles when gas density exceeds a "critical 
density" [see@]. Accretion of gas onto stars is calculated using the Bondi-Hoyle rate, which is 
capped at Eddington. Our units are M u = 3.5 x 10 6 M Q , the mass of Sgr A*, R u = 1.2 x 10 17 cm 
~ 0.04 pc, equal to 1" when viewed from the ~ 8 kpc distance to the GC, and t u = l/£l(R u ), 
the dynamical time evaluated at R u , approximately 60 years. We use R to signify distance in 
physical units and the dimensionless r = R/R u throughout the article inter-exchangeably. 

3. Initial conditions 

Our initial conditions all comprise a collision between two gas clouds at the edge of the inner 
parsec of the GC. The specific parameters for each run, labelled SI to S6, can be found in Table 
1. Each cloud is spherically symmetric and of uniform density, containing less than 1% of the 
mass of Sgr A*. The clouds are mainly composed of molecular hydrogen, using /i = 2.46. The 
initial temperature of the clouds is arbitrarily set to 20K. A more realistic model could have 
included a turbulent velocity and density field, but for practical purposes we limit ourselves to 
a smaller set of input parameters for this first study. 



Table 1. Initial conditions of our simulations. The meaning of the symbols in the Table are: (3 
is the cooling parameter, vi, v 2 are the initial velocity vectors of the two clouds, respectively; 
r pe and r ap are the pericentres and the apocentres of the two clouds; e is the eccentricity of their 
orbits; 9 is the angle between the orbital planes of the clouds, and b is the impact parameter. 
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We define the primary to be the larger cloud, with a radius of R\ = 0.2 pc and M\ = 
3.4 x 10 4 Mq, whilst the secondary has a radius of R 2 = 0.172 pc with M 2 = 2.6 x 10 4 M . The 
initial positions of the cloud centres are the same for all the simulations, and are r\ = (25, 0, 0) 
and r2 = (22,6, 7) in dimensionless units. The initial density of the clouds is slightly above the 
tidal shear density at their initial locations. 

In all the tests the absolute velocity of cloud 1 is of the order of the Keplerian velocity at 
that radius (and would be exactly circular in the absence of the cusp, M cusp = 0). Orbits in a 
cusped potential however are not Keplerian. The initial specific energy and angular momentum 
I of a particle can be used to find the pericenter, r pe , and the apocenter, r ap , of the orbit. The 
orbit's eccentricity e is then defined via 

\ + e _ r _2£_ ^2) 

1 e Tp e 

For any hyperbolic orbits, the e is capped at 1. The orbit of cloud 1 is slightly eccentric, with 
pericentre and apocentre of 25 and 31.5 respectively, and with e = 0.12. The initial orbit 
trajectory of the second cloud is varied between the tests to cover a small range of possibilities. 
The parameters for this trajectory in terms of pericentre, apocentre and eccentricity are given 
in Table [TJ The collision itself is highly supersonic. 

We run tests with cooling parameter f3 = 1 and /3 = 0.3 (see equation [1]) . These values are low 
enough so that fragmentation would occur if and when regions of the gas became self-gravitating 



13j;ll4|. Since the faster cooling runs were expected to require on average shorter timesteps the 



accretion radius for /3 = 0.3 was set to r acc = 0.33 whilst for j3 = 1 this was set to r acc = 0.06. 



4. Numerical simulations 

We shall now present some of the results of our simulations. In all of our runs, the clouds un- 
dergo an off-centre collision at time t ~ 10 (~ 600 years). As the cooling time is longer than the 
collision time, t co \\ ~ (Ri + i?2)/(|vi — V2I), the clouds heat up significantly and hence initially 
expand due to thermal pressure. This thermal expansion modifies velocities of the different 
parts of the clouds by giving gas thermal velocity "kicks" . The net result is a distribution of 
gas velocities that is much broader than we would get if the shock were isothermal. 

Gas dynamics 

The collision's product is not self-gravitating (until it has cooled down again) and so is easily 
sheared by the tidal field of Sgr A* . The collision and the resulting mixing of the clouds leads to 



angular momentum cancellation in shocks of some parts of the gas. Regions of gas that acquire 
smaller angular momentum infall to the respective circularisation radius on the local dynamical 
time. A small-scale disc around the black hole is thus formed on this timescale. Regions of 
the clouds that did not directly participate in the collision retain more of their initial angular 
momentum. These regions are sheared, then cool into filaments of length comparable to the 
initial sizes of the clouds' orbits. Parts of these filaments collide with each other or with the 
inner disc if the pericentres of their orbits are small enough. The inner disc therefore gains mass 
in an asymmetric manner, possibly resulting in strong warping and a change in its orientation 
with time. 

As the mass of the disc increases, it circularises, cools and undergoes gravitational collapse, 
with high density gas clumps being formed. Sink particles are introduced inside these clumps 
and allowed to grow in mass via gas accretion. We therefore see a significant amount of star 
formation, and in all our simulations, although the precise distributions differ, we always end 
up with at least two distinct stellar populations; in the disc and in the outer stream(s). 

T -7—^ j 



Figure 1. Gas eccentricity defined on 
radial shells function of the shell's 

radius at different times. Black, blue and 
red curves corresponding to time t = 50, 
100 and 250. Note that initially only gas 
on near-plunging high eccentricity orbits 
arrives in the innermost region, but with 
time gas circularises to e <C 1. 

Figure 1 shows the profiles of orbital eccentricity defined on radial shells plotted for several 
snapshots. These profiles show that soon after the collision the inner part of the computational 
domain is dominated by gas on plunging - high eccentricity - orbits. At later times, shocks force 
the gas to circularise at r ^ a few. Gas circularises faster at smaller radii as the orbital time is 
shortest there. As time goes on, nearly circular gaseous orbits are established in the innermost 
few arcseconds, whereas more eccentric ones dominate at larger radii. 

Star formation 

Stars in all of our runs form in both the disc and at least one filament. In the disc, stars form first 
in the inner arcsecond at t ~ 900 and later at a radius of ~ 5-8" at t ~ 1700. In the filament (s), 
star formation is approximately co-eval with the inner arsecond population and occurs at a 
radius of ~ 15-25". The disc and filament stellar orbits are inclined to each other by ~ 60°. 
Mass functions in the disc are top-heavy, whilst the filament (s) form entirely low- mass stars 
(0.1-1 Mq). However, as outlined in £}2j the mass spectrum of stars formed depends on several 
poorly constrained parameters of the simulations [e.g., see|l§]. Indeed, in the faster cooling runs, 
the stars in the outer filament(s) form mostly in a clustered mode. We expect that addition of 
radiative transfer and feedback would supress artificial fragmentation in these clusters and lead 
to a much more top-heavy mass function of the stars formed in the filament (s). 

All three stellar populations can be seen in Figure El We distinguish the locations at which 
they form by "inner disc" (within the inner 1-2"), "mid-range disc" (5-8") and "filament" (self- 
explanatory). The orbits in the inner disc population are almost circular, e ~ 0.05, whilst the 
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Figure 2. Gas surface density and stars for simulation SI at time t = 1955. Left: the whole 
of the central parsec, showing the inner disc, mid-range disc, and filament stellar populations. 
Right: a zoomed-in plot of the inner disc. Sgr A* is located at (0,0), and the line of sight is 
along the y direction. Stars are colour-coded by mass: green (0.1-1 M Q ), cyan (1-10 M Q ), dark 
blue (10-150 M ), magenta (> 150 M ). 



mid-range disc population has an eccentricity of e ~ 0.2. Both of these are therefore in good 
agreement with the clockwise feature of the observations [ij. The orbit of the filament popu- 
lation is also eccentric, with e ^ 0.2. This is in reasonable agreement with the corresponding 
counter-clockwise feature of the observations, although the eccentricity of the latter is considered 
to be significantly higher with e ~ 0.8 0]. 

Infall of gas onto inner disc 

When the feeding into the inner 1-2" is reasonably steady-state, and cooling is long, stars form 
here once the gaseous inner disc has settled into a stable orientation. The resulting stellar 
distribution of is that of a geometrically thin disc [H/R < 0.1), a distribution that is broadly 
consistent with the observed orbits of young massive stars in the clockwise disc [§]. 

In some of our runs, however, the inner disc undergoes a complex evolution in angular mo- 
mentum space. When the feeding of the disc due to infalling gas is intermittent in its distribution 
of angular momentum e.g. due to a more grazing initial collision, the newly arriving gas tweaks 
the disc orientation significantly. However, if stars are introduced in the disc whilst its midplane 
orientation is still changing (as is the case for the (3 = 0.3 runs), these stars remember the 
"old" orientation of the gaseous disc in which they were born whilst the gaseous disc evolves 
to quite a different orientation. Possibility of this effect taking place was suggested based on 
analytical arguments by 0]. These authors found a "critical rotation timescale", below which 
the disc will leave the stars behind as its midplane rotates. Indeed, in one of our tests with 
fast cooling and a grazing initial collision, the inner disc rotates on a timescale much shorter 
than this critical value, and this coupled with several star forming events in gaseous discs of 
different orientations leads to a geometrically thick stellar disc (H/R ~ 1) in the inner arcsecond. 
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Figure 3. Stellar surface density, 
£*, versus projected radius, for 
the run SI at t= 1955. Different 
colours show the three different 
orientations, along the z, x and 
y axes (black, red, and blue 
respectively) . 



Radial distribution of stars 

Interestingly, the most robust result from all of the 6 runs completed is the fact that distribution 
of stellar mass in an annulus, dM* = Y^ if (R)2TrRdR, versus radius roughly follows the law 

oc -L . (3) 

We show these distributions for run SI in Figure HI To facilitate comparison with observational 
data, we use projected radius p instead of the proper 3D radius r, which is not directly known in 
observations of young stars in the GC [e.g., |2|. As the viewing angle of the stellar system mod- 
elled here is arbitrary, we chose to plot the distributions along the three axes of the simulations, 
so that p = \J x 1 + y 2 (black symbols), p = y 7 y 2 + z 2 (red), and p = \J x 2 + z 2 (blue). The lines 
in the Figure show R~ 2 (dotted) and i?~ 3 / 5 (dashed) power laws. The latter corresponds to 
that predicted by a not self- gravitating standard accretion disc theory. This result is remarkable 
given very large differences in the 3D arrangments of the stellar structures in our runs. The 
1/R 2 law is in a good agreement with the observed distribution [3]. 

Sgr A* feeding during and after the cloud collision 

Figured] shows the mass accretion rate for the central black hole for simulation SI. Formally, our 
simulations provide a sustained super-Eddington accretion rate ( ^ 0.03 M yr _1 ), in some cases 
for over 10 4 yrs, and at a fraction of that for up to 10 5 years. However, we do not resolve gas 
dynamics inside the accretion radius, where the material should form a disc in which accretion 
proceeds viscously. The viscous time scale, t v i sc , depends on the temperature in the disc midplane 
and the viscosity parameter, a [3]. The midplane temperature in the inner arcsecond is around 
lO 3 ^ in both the standard [non self-gravitating, e.g., 1H and the self-gravitating regimes (lH |. 
yielding H/R ~ 0.01: t v isc = 6 x 10 years \(R/W0H) 2 r^ 2 , where a .i = a/10 and r is m 
our code units (one arcsecond). With these fiducial numbers, the viscous time scale coincides 
with the age of the young massive stars in the GC [e.g., 0; If i v isc *C a few million years 
(e.g., if a = 0.1, t v ; sc ~ 10 6 years at 0.3"), we expect that gas would have mainly accreted onto 
Sgr A* by now. This accretion rate would be a significant fraction of the Eddington accretion 
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Figure 4. The time evolution of the 
accreted mass onto the central black 
hole (solid line) and the corresponding 
accretion rate (dotted line) compared 
to the Eddington limit for Sgr A* 
(dotted red line), for simulation SI. 



rate. Standard disc accretion would then generate as much as 2 x 10 erg M acCi 3 of radiative 
energy, where M acCj 3 is the gas mass accreted by Sgr A* in units of 10 3 Mq. A similar amount of 
energy could have been released as energetic outflows. There is currently no evidence for such 
a bright and relatively recent accretion activity of Sgr A*. 

In the opposite limit, i.e. if a -C 0.1, and i v i sc a few million years, the gaseous disc should 
still be there as self-consistent modelling predicts that stars should not be forming within region 
r ^ 0.3" — 1" [16]. This would contradict the data, as no gaseous disc is observed now [3; 111]. 
Therefore, our simulations seem to over-predict the amount of gaseous material deposited in the 
inner ~ 0.3". We take this as an indication that a better model would perhaps involve a gaseous 
cloud of a larger geometric size, thus shifting all spatial scales outwards. 



5. Discussion 

In approximate order of robustness, our key results are as follows: 

(i) Formation of a gaseous nearly circular disc in the inner region of the computational domain 
is common to all runs, as is the ensuing formation of stars on similar circular orbits. This 
is natural as the dynamical time in the innermost disc is only ~ 60r 3 / 2 years. Conversely, 
the outer gaseous stream becomes self-gravitating much faster than it could circularise, 
and hence orbits of stars in that region are more eccentric, in reasonable agreement with 
observations 0]. 

(ii) Radial distribution of stellar mass closely follows the observed oc 1/R 2 profile of the disc 
stellar populations 0]. This was observed for all the runs, although we expect these results 
to change if colliding clouds moved in similar directions, significantly reducing the angular 
momentum cancellation in shocks and the thermal "kick" velocity due to the shock. 

(iii) Runs with a comparatively long cooling time parameter, [3 = 1, lead to kinematically less 
dispersed stellar populations than runs with faster cooling. As a result, the longer the 
cooling time, the more closely the resulting stellar system can be fit by planar systems in 
velocity space. The innermost stellar disc is then reminiscent of the observed clockwise thin 
stellar system [3]. Rapidly cooling runs produce clumpier gas flows that lead to significant 
gaseous disc orientation changes, producing less coherent discs; such geometrically thick 
systems are incompatible with the observations. 



(iv) With the chosen initial conditions, faster cooling promotes survival of the gaseous streams 
corresponding to the orbits of the original clouds. These streams fragment and form stars 
mainly in a clustered mode, although this is expected to depend on the details of the 
radiation feedback from young stars, which is not modelled in this set of simulations. 

The observed well defined, flat, geometrically thin and almost circular clockwise stellar system 
0] is best created via a gentle accumulation of gas. Several independent major gas deposition 
events lead to a too warped disc, and/or mixed up systems consisting of several stellar rings or 
discs co-existing at the same radius. To avoid this happening, the inner disc must be created on 
a timescale longer than the critical rotation time, which is estimated at t cr ~ few xlO 4 years. 

Deposition of gas onto the inner disc will most likely take place over the collision time, 
t C oii ~ R c \/v c \, where R c \ and v c \ are the cloud's size and velocity magnitude. We therefore 
require that the collision itself be more prolonged than t CT . Estimating the velocity of the cloud 
at v c \ ~ 150 km/sec, which is of the order of circular velocities in the GC outside the inner 
parsec, we find t co u = R c \/v c \ ~ 10 4 years -R c i,pc 5 where -R c i iPC is the size of the cloud in parsecs. 
We hence require the cloud to be larger than a few parsecs to satisfy t co \\ ^ t cr . Note that this 
size is not necessarily the original size of the cloud if the cloud gets tidally disrupted before it 
makes the impact. In the latter case we can take R c \ to be the radial distance to the centre of 
the Galaxy at which the tidal disruption took place. Finally, the location of the collision should 
not be too far from the central parsec, or else too much angular momentum would have to be 
lost to deposit a significant amount of gas in the ~ 0.1 pc region. 

In addition, the radial distribution of gas and stars in our simulations is too compact, 
contradicting the observations (no He-I stars inside the inner arcsecond and no gaseous disc 
there either). We thus believe that a realistic scenario would be a GMC of the order of a few 
parsecs in size striking the CND a few parsecs away from Sgr A*. Alternatively such a cloud 
could self-collide if the impact parameter with respect to Sgr A* is small enough, but the cloud 
would need to be very structured, e.g., essentially consist of several smaller clouds or filaments. 
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